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Abstract 



O ■ We discuss a model in which the fundamental scale of gravity is restricted to 

^ ! 10~^ eV. An observable modification of gravity occurs simultaneously at the Hubble 

distance and at around 0.1 mm. These predictions can be tested both by the table- 
top experiments and by cosmological measurements. The model is formulated as a 
' brane-world theory embedded in a space with two or more infinite-volume extra di- 

mensions. Gravity on the brane reproduces the four- dimensional laws at observable 
^ ! distances but turns to the high-dimensional behavior at larger scales. To determine 

^ I the crossover distance we smooth out the singularities in the Green's functions by 

taking into account softening of the graviton propagator due to the high-dimensional 
operators that are suppressed by the fundamental scale. We find that irrespective 
of the precise nature of microscopic gravity the ultraviolet and infrared scales of 
gravity-modification are rigidly correlated. This fixes the fundamental scale of grav- 
ity at 10~^ eV. The result persists for nonzero thickness branes. 



1 Introduction 



The law of gravitational interactions is tested experimentaly at distances above 0.1 
millimeter [|l|. Hence, gravity could change its behavior at a scale which is as low 
as: 

Mf P ~ (0.1 mm)-^ ~ 10"^ eV . (1) 

A theory of gravitation which predicts a deviation from the conventional law above 
M*, such that gravity is still negligible above this scale compared to the gauge 
interactions, cannot be ruled out at present. 

Another important piese of data which is relevant for gravity comes from recent 
astrophysical measurements. The data indicate that the expansion of the univesre 
might be accelerating The existence of the cosmic acceleration raises yet an- 
other interesting question on possible infrared (IR) modification of gravity. Indeed, 
an alternative view might be that the observed acceleration of the universe is not 
because of the vacuum energy, but due to the modification of the laws of gravity 
at the Hubble distance. Therefore, there might be another scale in the theory of 
gravitation 

H^^ ~ 10^^ mm , (2) 

where Hq denotes the present-day value of the Hubble parameter. 

The aim of the present work is to discuss a model which predicts the modification 
of gravity at distances of order 10^^ mm and as a consequence of this, the modifica- 
tion at 0.1 mm ensues. As we will see, these two scales are intimately related. As a 
result, the model restricts the value of the fundamental scale of gravity to 10^^ eV. 

A brane-world model, in which extra dimensions have an infinite- volume and the 
4D laws of gravity are reproduced on a brane due to the induced Einstein-Hilbert 
term, predicts the modification of gravity at very large cosmological distances [Q. 
The 4D Einstein-Hilbert term "shields" a brane observer from strong bulk gravity . 
Consequently, the observer detects the conventional 4D gravity all the way up to 
large cosmological distances. At those scales the effects of infinite-volume extra 
dimensions take over and the behavior of gravity changes. In a 5D model one finds 
the following crossover scale [§: 

where Mpi is the 4D Planck mass and M^, denotes the fundamental gravitational scale 
of a 5D bulk theory. For large distances, r 3> r[!^^^\ gravity is five-dimensional, 
while it is four-dimensional when r -C r^^^^\ The laws of the 4D gravity are valid 
all the way down to the distances of order However, below this length scale 

the effective theory of gravity breaks down. As a result, there is a lower bound on 
the scale which comes from the accelerator, astroparticle, and cosmological data, 
that is > 10-3 eV M. 
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Therefore, in a 5D theory one finds the laws of 4D gravity only in the interval: 



±<r<^. (4) 

This scenario leads to novel astrophysical and cosmological consequences and predic- 
tions which were found in Refs. |5Hl^. In fact, it allows for the accelerated universe 
with no cosmological constant which is consistent with the data Some of these 
results will be discussed below. A similar phenomenon for photons was studied in 
Ref. ||15| (see also Ref. In these models the behavior of gravity and gauge fields 

are alike. 

Note that in order for the large distance modification of gravity to have relevance 
to present day cosmology and astrophysics the crossover scale should be at, or below, 
the observable size of the universe, r^^^^^ < 10^^ mm. This imposes the following 
restriction > 10 MeV. Therefore, in this scenario, gravity will simultaneously 
be modified at the Hubble distance and at (10 MeV)^^ ~ 20 Fermi. 

In the present paper we will argue that the higher dimensional {D > 6) gen- 
eralization gives rise to a different dependence of the crossover scale on Mpi 
and M*. In particular, for the modification of gravity at around 10^^ mm we get 
~ 10~^eV, and vise versa. The fundamental scale of gravity is bound to be 
10^^ eV. Thus, this model makes very restrictive predictions on the modification 
of gravity which are testable in table-top experiments and simultaneously can be 
testable in astrophysical and the cosmological measurements. 

There are additional reasons why this model can be interesting to studyQ The 
scale of 10~^ eV is very close to the scale of the vacuum energy density 0. Although, 
this similarity might be a coincidence, nevertheless, it is rather intriguing to argue, 
following Ref. |1^ , that the scale of the vacuum energy density in the Universe might 



be related to the ultraviolet cutoff of the gravitational theory. Moreover, theories 
with infinite-volume extra dimensions might shed new light on the cosmological 
constant problem since the bulk SUSY can be preserved in these models unbroken 
2^,^, p (for discussions of these models in string theory see Refs. [p^-p^). 



The generalization to the D = {A + A^) -dimensional models with A^ > 2 is not 
completely trivial, however. Certain additional singularities emerge in this case 
and they should be dealt with some care. The presence of these singularities is 
related to the fact that the Green's function of the d'Alambertian in the space 
transverse to a 3-brane blows up at the origin for N > 2, while it is regular for 
A^ = 1. One of the purposes of the present paper is to deal carefully with these 
singularities. 

In what follows we will cast the singularities into the following two classes. 

• Type I: These are the singularities which are related to the fact that the brane 
is approximated by a delta-function type source. They are present as in the 
A^ = 1, as well as in A^ > 2 cases. 



^For models with warped geometries and modification of gravity at large scales see Refs. ||l7| , |l8| . 
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• Type II: These are the singularities which are present for the N > 2 case; they 
are related to the fact that the transverse Green's function of the d'Alambertian 
in iV > 2 dimensional space blows up at the origin. 

Although these two types of singularities have a common origin, nevertheless, it 
is convenient to discuss them separately in order to compare the N = 1 and N > 2 
models. 

Let us start with the Type I singularities, i.e., those associated with a vanishing 
brane width. Suppose the brane had a finite thickness (we call it A)0. What would 
be the effect of this width on a low energy observer in the 4D worldvolume theory? 
One effect is the presence of a massive state (or states) which correspond to the 
fluctuation of the brane transverse width (the "breathing mode"). The mass of this 
state is of the order of the inverse brane width. This state, being localized and 
gravitationally interacting, will certainly renormalize the graviton kinetic term on 
the 4D worldvolume, so that the generic value of the coefficient in front of the 4D 
graviton induced term is of the order A~^. Other than that, the presence of such 
a heavy state is irrelevant for low-energy 4D physic^. Therefore, from the point of 
view of a low-energy observer this state can be "integrated out" . This is equivalent 
of "integrating out" the brane width. Hence, in the low-energy theory the brane 
will look as a delta-function source. Since we expect that the physics is smooth and 
continuous, the calculations at low energies, i.e., at £^ -C A~^, in a theory with the 
delta- function brane should produce the same result (up to the 0{EA) corrections) 
as calculations in a theory where the ultraviolet resolution of the brane is manifest 
due to a nonzero brane width. This is indeed what we will find below. 

Let us now turn to the Type II singularities. They will be taken care of separately 
from the singularities associated with the delta-function brane. In fact, as a first 
step we could keep the brane as a delta-function source and try to remove the Type 
II singularities. For this one could introduce the "e shift regularization" as it was 
done in 0, or one could use the rigid cutoff in the bulk as in Ref. However, 
a more justified way (that automatically preserves reparametrisation invariance) 
to deal with this problem is to introduce higher dimensional HD operators as was 
proposed in Refs. p2[ and [Q. From the physical standpoint, the HD's are remnants 



of UV physics which is generically unknown. Thus, it is natural that they provide 
regularization of certain UV singularities in the low-energy theory. This is the way 
which we will follow in the present paper. 

The paper is organized as follow. In section 2 we study the induced terms on 
a thick brane. We will show that these terms in general are nonlocal and have no 
factorizable form in terms of the worldvolume and transverse coordinates. These 
nonlocal expressions, however, can be expanded in an infinite series of local terms. 



^ Below in sections 3-6 we will be dealing with a brane the transverse width of which is much 
smaller than the fundamental distance scale in the bulk. Such an unusual brane can exist even in 
4D theory, see discussions below and in section 3 of Ref. [Q . 

■^It also renormalizes the brane tension. 
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The expansion parameter is the width of the brane. The leading term coincides with 
the 4D Einstein-Hilbert term. In section 3 we setup a scalar field theory model which 
mimics the properties of a theory with the induced kinetic term on a thick brane. 
At this point we fine-tune the brane tension and the 4D cosmological constant 
to make the calculation analytically treatable. In section 4 we use this model to 
calculate the two-point Green's function. First we consider the delta-function brane 
but introduce the HD operator which smoothes out the Type II singularities in the 
expression for the Green's function. We show that the behavior of the gravity is 
qualitatively similar to that of the 5D theory: at short distances the theory behaves 
as four-dimensional and at large distances it approaches the (4 + A^)-dimensional 
regime. The expression for the crossover distance in this case is different from the 
one found for the 5D theory. In section 5 we consider the case when the brane width 
is kept fixed and calculate the corresponding Green's functions. We show that the 
qualitative behavior described above holds unchanged. The limiting transition to 
the delta-function brane is smooth. There are no ultra-local interactions lEBI in 



the bulk and no new physical scales p2| arise in these calculations. In section 6 
we discuss the phenomenological constraints on the crossover scale beyond which 
gravity can change its regime. In section 7 we consider a brane with a nonzero 
tension and give an estimate for the crossover scale in this case. Conclusions are 
presented in section 8. 



2 Induced Terms on Smooth Branes 

In this section we study how the induced terms arise on a smooth solitonic brane. We 
will argue that if a nonzero brane width is kept, the induced terms are nonlocal and 
not factorizable with respect to worldvolume and transverse dimensions. However, 
in the low-energy approximation they can be expanded in a power series of local 
terms with the leading term coinciding with the 4D Einstein-Hilbert action. 

Our consideration in this section is rather general and applies to any brane- 
world model (with a non-conformal-invariant worldvolume theory). Let us start 
with a (4 + A^)-dimensional action which contains a graviton Gab and a scalar 
field H which makes the brane, and other fields \1/ which do not participate in the 
formation of the classical background^: 

S = J d^x d^y \fG H, ^) . (5) 

Here £, denotes the total Lagrangian density which is a function of the fields and 
their derivatives which are suppressed in (^. The bulk graviton is decomposed as 
follows: 

GAB{x^,yi), A,B = Q,l,. + N- /x, i/ = 0, 1, 2, 3; i, j = 4, 5, .., 3 + AT . (6) 

^ All the derivations in this section can be generalized straightforwardly for more complicated 
classical backgrounds. 
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We study the Lagrangian which yields a classical solution in a form of a three- 
brane. This solution depends on the coordinates i/i but does not depend on x^. 
Many examples are known in various dimensions. Let us split the graviton and 
scalar fields in their classical parts and fluctuations: 

Gab = G%{y) + HAB{x,y) , G%{y) ^ A\y) r^,, , G% = , (7) 

U{x„y,) = n^\y) + aix,) f{y) . (8) 

Where in (|^) we use the parametrisation of the {/uz/} component of the background 
metric in terms of the function A'^{y) which is typical for the branes that preserve 
the 4D Poincare invariance on the worldvolume. 

The scalar field 11 has modes which are localized on a brane. The most obvious 
ones are the Goldstone bosons associated with spontaneously broken translation 
invariance in the yi directions. These particles are massless and are derivatively 
coupled to matter fields on a brane. For these reasons they are not relevant for 
our discussions below. In addition, generically there are massive localized modes on 
a brane. For simplicity we consider below only a single mode. As we mentioned 
before, the latter corresponds to the fluctuations of the transverse width of the 
brane. Thus, its mass is proportional to the inverse brane width. Let us discuss 
this mode in detail. Since it is localized on a brane there should exist an effective 
four-dimensional Lagrangian for it. To derive this Lagrangian let us start with the 
action for 11: 

J d^x d^'y VG dAlldBll - V{U)^ . (9) 

Here is a potential which is responsible for the existence of the brane. Let us now 
substitute and (||) into @ and keep truck of quadratic fluctuations (interactions 
will be discussed later). In this approximation we find: 

(The prime denotes differentiation w.r.t. 11. The 4D indices here and below are 
contracted by the fiat space metric rj^y and those of the extra coordinates by Gf^) 
The quadratic term in a is a mass term. The value of the mass depends on the 
exact expression for /. For instance, by definition of a localized Goldstone particle 
its profile / is such that the mass term in ([T0|) vanishes. On the other hand, the 
profile / for a breathing mode is different, it satisfies the equation |: 



fd,[ ^G^\y)d^f) _ p 



f V"(ff') + L = _J . (11) 



^Note that functions / for different modes are orthogonal to each other, that is why we can 
treat them separately. 
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Here the parameter M is proportional to the inverse brane width: 



M ~ 4 • (12) 



Using ([TlD we obtain the following expression for the quadratic part of the low 



energy 4D action of a breathing mode: 

We see that the y dependence can be integrated out. This is because for large y the 
localization function for local defects is exponentially decreasing, (~ exp(— |?/|/A)). 
This has to make the y integral in (|13D to converge for a breathing mode even though 
the integrals 



^lG^)d''y and / ^fcF{y)/A\y)d''y 



might not be convergent (the latter case corresponds to the brane models where 
there is no localized graviton zero mode). 

Let us now turn to the gravitational part of the action. We will consider a 
general case when the graviton is not necessarily localized. If so, the bulk gravity 
action cannot be integrated w.r.t. y to obtain 4D kinetic term for 4D gravitons. 
Nevertheless, the 4D laws of gravity can be obtained on a brane due to the induced 
terms [0,^ to the discussion of which we turn now. 

Let us look at the interaction of bulk gravity with the localized field a. We 
concentrate first on the {^iv} components of the interactions since these are the 
ones that determine the 4D induced terms relevant to us. Using (0),(|D and (||) we 
derive the interaction Lagrangian for the {/iz^} part: 



j d'x d'^y ^fG^)^^r'V^^H^p{x.y) X 



A\y) 

]^dxa{x)d^a{x) - ^^'M^jj . (14) 

The quantity in second line in (|14D is nothing but the energy momentum tensor for 
a with the mass M. From this we see that the interaction of a bulk graviton with 
the localized mode can be rewritten in a purely 4D form. Indeed, rescaling the field 
0" — > cr/(/ d^y^ G^^{y) P{y) / A^iy))^^"^ we get the canonically normalized 4D kinetic 
term for the field a in (|13]) and its interaction with the bulk gravity, according to 
Eq. ([l^ , takes the following form: 

5int = / d^'x h^^{x) T;,(x) , (15) 
where the field h is defined as follows: 



, ^ d^y[^G^) f\y)H,.{x,y)/A\y)] 

hf,u{x) = . . (16) 

/ d^y [.Jg^) P{y)/A^{y)] 
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Thus, the interaction Lagrangian of higher dimensional gravity with a locahzed 
mode can be given a purely 4D form (^) in spite of the fact that the bulk graviton 
might not be localized on a brane (i.e., the bulk graviton kinetic term might not be 
integrable w.r.t. y). In what follows we will study loop effects due to this interaction. 

Before we turn to these issues we would like to digress for a moment and make 
two comments. So far we discussed only a single localized field. Suppose now there 
are n localized fields ak ,k = 1,2. ..,n, with the localization functions fk which in 
general can be different. Then, the interaction vertex of each of this state with the 
bulk gravity can be presented in the form of Eq. (|T5|) . On the other hand, since /'s 
are different and the expressions for /'s determine h^^{x) in (p!6[), then each of these 
localized states will interact with different "effective 4D gravitons" (p^Gf) . However, 
the effects that discriminate between these gravitons are suppressed by the brane 
width and are negligible as we will see below. For this let us perform the expansion: 

H^,{x,y) = H^,ix,0) + y'diH^,{x,0) + ^ y' y^ d,d,H^,{x,0) + ... . (17) 

We substitute this series into (|16]) and then in (plSj) for each localized field. As a 
result, for each field we get infinite number of terms first two of which are: 

^'^'^(^,0)T-(x) / \/G%)^d^, (18) 

m^^{x,o)T^^;;ix) Jy' \/G^)§^^d''y^ (19) 

and so on, with increasing powers of derivatives. The first term (^) gives the 
universal interaction for all the different fields. This can be seen by performing the 
rescaling of the fields to bring their kinetic terms in ( pTSf ) to the canonical form. Then, 
the couphngs of all the localized fields to H^^{x, 0) are identical. In fact, the result 
of this universal term is what we obtain in the delta-function limit for the brane, i.e., 
when p 5^^\y). This is the dominant contribution. The higher terms, on the 
other hand, are not like this. Each higher term (including the one in ([T9|)) can be 
thought of as a new field didj...dk H^u{x, 0) from the 4D point of view. However, as 
we can see from (0), the couplings of these fields to the energy-momentum tensor 
contain extra powers of y under the integral. Since the effective region of integration 
is determined by the "width" of the function /, then each additional power of y will 
translate into an additional power of A. Therefore these couplings are suppressed 
by extra powers of the brane width A. For instance, the term in (p!9| ) can give rise 
to the contributions containing extrinsic curvature and its derivatives. However, 
these terms are suppressed compared to the leading ones by the powers of the brane 
width0. 

^In Ref. p3|| it was found that on a D-brane in bosonic string perturbation theory the induced 
extrinsic curvature terms are of the same order as the Einstein-Hilbert term and both are suppressed 
by powers of string couphng constant. However, the nonperturbative effects, e.g., those due to 
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Therefore, we conclude that in the low-energy theory the operators which would 
violate universality are suppressed by powers of A. When the brane width is small 
enough, of the order of 1/Mpi, these effects are unobservable at low energies. 

After this digression we return back to the case of one localized scalar field and 
try to study it as far as possible. 

Having the interaction vertices derived let us look at the loop diagrams where 
two external graviton lines are attached to the loops in which the localized fields are 
running. The vertices in these diagrams are described by (p!5D. One can think of 
this interaction as a purely 4D interaction of Tl^^{x) with the effective 4D graviton 
h^y{x). The loop diagrams with two external h^i, lines give rise to the following 
term in the effective action on a brane: 

KM d^''^^{p,M)Kp{-p) , (20) 

where 0^'^°'^{p, M) is some function of the external momentum p and the particle 
mass M (in fact, the UV cutoff in this case is M ~ A~^). In the one loop approx- 
imation this function can be expressed in terms of hypergeometric functions (see 
e.g, [0). This leads to a nonlocal interactions in the effective Lagrangian in the 
coordinate space. The explicit form of the nonlocal coordinate-space Lagrangian is 
hard to present. However, at the momenta p <^ M we can perform an expansion of 
the form-factor 0^^°"^ in powers of p^ /M"^ . The leading term takes the form: 

hi,v{p) (vi^aVuf^p'^ - VfiuVa/Sp"^ " Vt^aPuP/B 
-VupP,iPa + VapPfiPu + Vt^iyPaP/s) Kp{-p) ■ (21) 

This is a quadratic approximation to a reparametrisation invariant graviton kinetic 
term in a 4D low-energy effective theory of gravity (for treatment of general relativity 
as an effective field theory, see Ref. p8| ). 

Let us now turn to the {ij} components. These will induce the brane-kinetic and 
mass terms for the components Gij. The latter, look as scalars ( "graviscalars" ) from 
the point of view of a braneworld observer. However, the mass terms (potentials) 
of these states on the brane might not be protected by any symmetries and acquire 
the mass of the order of the UV cutoff of the worldvolume theory. Such states 
cannot mediate forces which would compete with gravity at observable distances. 
In addition, even if they stay massless, they should decouple from the 4D matter at 
short (observable) distances in analogy with the phenomenon found in Ref. . 

Finally, we turn to the {i/x} components. From a 4D perspective these are 
graviphotons. In the linearized approximation they are derivatively coupled to 
conserved currents. These terms cannot give rise to gravity competing forces at 
observable distances and are irrelevant for our considerations. 

the inverse brane width studied above, are expected to give the dominant contribution to the 
worldvolume Einstein-Hilbert term. The terms with extrinsic curvature will in general be induced 
on a brane as well. It is interesting to study the effects of these terms. 
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3 A Scalar Field Model 



To summarize the results of the previous section we write down the gravitational 
part of the action which includes the lowest dimensional derivative terms in the bulk 
and those induced on the brane: 

S = Ml+^ J d^x d^y ^\ n{G) + Ml^ J d^x ^\ R{g) + other terms . (22) 

Here, the induced 4D term should be understood as an expansion in perturbations 
defined as follows: 

gfiu{x) = 7]^u + hf,y{x) , (23) 

where h^y{x) is related to the bulk graviton G via (p!6l)PT. The first term in (^) is 
the bulk kinetic term for a graviton and the second term is a brane-induced kinetic 
term for the very same graviton. 

In the previous section we found that in general the induced terms are nonlocal. 
The term which is given in (^) is just a first term in the expansion of that nonlocal 
expression in powers of the brane width A. However, as long as we discriminate 
between h^y{x) and H^u{x^ 0), the induced term in ( p2|) itself contains the effects of 
the same order, i.e., those which are proportional to A. This means that for any 
nonzero A the action ( ^2]) is not complete. The additional terms which vanish in 
the A — > limit are missing in (^2]). Nevertheless, the latter action can be used to 
check the limiting transition from the finite thickness brane to the delta function 
brane since in the limit A — the missing terms do not contribute. Therefore, we 
will be able to check whether the action (^) gives rise to the results which in the 
limit A — > turn smoothly to the ones obtained with a delta function brane []3|,|5|,p!5 
answering some question raised in Refs. P^,pB|. 



Our goal is to study the interactions between sources that are located on the 
brane and which exchange gravitons described by the action (^2]). The value of the 
induced constant is crucial. In order to reproduce the correct 4D gravity on 
a brane this has to equal to the 4D Planck mass, Mf^^ = Mp^ |^,^. The question 
is how does one get such a big scale on a brane when the bulks scale M^, is rather 
low. There are some known possibilities for this which we will mention briefly. 
First, there exist certain branes, even in 4D theories, the width of which is much 
smaller than the fundamental length scale in the bulk (see discussions on this in 
Sect. 3 of Ref. |^.). In this case, the induced constant will have the magnitude 
much bigger then the bulk scale M^. More attractive possibility is to consider a 
brane worldvolume particle theory that has very high UV scale (such as the GUT 
scale) . States of the particle theory renormalize the gravitational constant and make 
it of the order of the UV cutoff. In general, the induced scale is determined by a 
two-point correlation function of the energy-momentum tensor of the worldvolume 
particle theory which is localized on the brane. This was discussed in detail in 



^Note that the induced metric on the brane is actually A^(A)r7^^. The constant factor A'^{A) 
rescales M^^; it is dropped for simplicity above. It will be taken into account in section 7. 
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Refs. 1^ (see also 0). As a result, the 4D gravitational Planck scale is a derived 
quantity. Thus, the relation Mjnd = Mpi is a definition of the 4D Planck mass Mpj. 
The next step is to calculate the two-point Green's functions for the gravitational 



action (22). This calculation can be straightforwardly performed for gravitons along 
the lines of [^. However, to give more clear derivation it is convenient to make the 
following simplification. In what follows we will suppress the tensorial structure since 
the latter does not change the issue of singularities. Although, there will be separate 
issues due to the tensorial structure (such as the ones associated with perturbative 



discontinuity and nonperturbative continuity in l/r^ advocated in Ref. |12]). Since 
physics of higher codimension cases is very similar to that of the 5D case, we expect 
that these issues can be addressed is a manner similar to what was found for the 5D 



case in Ref. |]T2[ based on earlier work (see also recent work |T^). Having this 
said, in what follows we substitute the graviton field by a scalar field $ and study 
the Green's functions for the latter. This allows to determine the force-law. Hence, 
the scalar $ should be regarded as a counterpart of a graviton. 

Until section 7 we will work with a brane which is placed at an orbifold fixed 
point and the tension of which is tuned to zero by means of the 4D cosmological 
constant. This can be applicable, for instance, to a non-BPS system of D-branes 
and Orientifold planes placed on top of each other at certain fixed points. The 
generalization to the nonzero tension branes will be seen in section 7. Before that, 
we should put = 1, A{y) = 1 and / P{y)d^y = 1. 

The action for the scalar field which mimics the graviton with the induced kinetic 
term on a brane in a space with extra dimensions is as follows: 



S = I d^xd^y M^+N 



+ Ml, J d'x dV dV f'(y') d,^x,j') d^<l>{x,j") . (24) 

Note that we have normalized the scalar field in such a way that it is dimensionless. 
The HD term, which is proportional to e, is included for the regularization of the 
propagator (see below). The parameter e is an arbitrary 0{1) constant. Let us 
mention here that the results obtained in the next sections are not bound to the 
form of the HD operator which we use. Any HD operator which smoothes out the 
UV effects will reproduce the same qualitative results. 

Below we discuss two separate cases. First we approximate the brane by delta- 
function and show that singularities are removed by the HD terms. In this case: 

f{y) = 5(^)(y) . (25) 

As a next step we consider a smooth nonzero thickness brane and study the tree-level 
propagator. For this we use the simplest ansatz: 

/^(^)- {oN-lyl>'A^: '''' 

where in the limit A ^ 0, ^ cxd and a^A^ = const. 
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4 A Delta Function Brane 



In this section we restrict ourselves to the consideration of a deha-function brane. 
The singularities will be removed by means of the HD operators. We adopt the 
expression ( p5D for the localization function and use the method of Ref. 0. The 
equation for the Green's function for a source located on the brane takes the form 



M, 



2+7V 



+ G{x,y; 0,0) 

= d^^^x) 6^''\y) . (27) 



To find a solution of Eq. ( pTj) let us turn to the Fourier-transformed quantities with 
respect to the worldvolume four-coordinates x^: 



Gix,y; 0,0) 



e^P" G (p, y) . 

(27r)' 



Then, in the Euclidean worldvolume space Eq. (p7|) takes the form 



M, 



2+N 



+ M^,6^''\y)p'jG{p,y) 

= s^'^Hy) , (28) 



where p = ^Jpf + pi + pI + pl denotes the magnitude of the Euclidean momentum. 
We are looking for the solution in the following form 

G{p,y) = D{p,y) B{p), 

where D{p,y) is defined as follows 

D{p,y) = 5W(y). 



(29) 



B{p) is some function which should be determined. Using the above expression one 
finds §: 

D{p,y) 



G{p,y) 
where D{p,y) is given by 

D{p,y) = 



Ml^p^D{p,0) + M2+^ 
d^q exp(zgy) 



(30) 



(27r)'' p^ + + t {p^ + qY/M^ 

Note that for > 2 the e/M^ term smoothes out the UV divergences in the 
expression for D{p,y). Indeed, without this term the expression for D{p,y) blows 
up at ?/ = 0. This divergence is not related to the presence of the delta function 
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brane as it is clear from Eq. (^). Below we keep the HD term (in the = 1 case 
the expression for D{p, 0) is finite even without the HD term P,|5|,p!5[].) 

As an example of a theory with higher than one codimensions we consider the 
case of = 3. The corresponding calculation can be done exactly. The results for 
N > 2 are similar 

For A^ = 3 the function D{p,y) is given by 

1 . . I 

D{p,y) = {exp{-p\y\)- exp{-\y\p)} ; p = \ p^ + — . (31) 

4vr \y\ V e 

At y = 0, we find 

which is finite. Note that, an action with HD terms has a consistent interpretation 
only if it is regarded as an infinite series in derivatives. Any truncation to a finite 
order can give rise to ghosts. The HD term should be treated as a {p^ / Ml) correction 
in the expansion. Hence, the above expression for D{p,y) makes sense only as an 
expansion in [p^ / Ml) in which only the first correction is to be kept. Combining 
this expression with Eq. (|30D , we get 



G(p,0) ^ ^ 



Mt^e^T^il +p[^Y] +Ml,p 



When the M^^p^ term in the denominator dominates this gives the 4D Green's 
function. The potential mediated by the scalar field $ on the 4D worldvolume of 
the brane is determined as 

V{r)= f G{t,x',y = 0;0,0,0)dt, 



where r = y xf + X2 + x^. Using Eq. (0), we get 

V{r) = r dp 2psin(pr) 1 ^33^ 



This potential behaves as that of a 4D theory when y/e/M^ -C r -C Mpi/M^ and 
behaves as the potential of a (4 + A^) -dimensional theory for bigger values of r. 
The similar calculation for the crossover distance can be done for other A^'s. The 
qualitative result is the same, i.e., the crossover distance is approximated by: 



(D>6) _ ^ 

Mr 



®Note that for > 3, in order to regularize the function D, we should include higher order 
operators beyond the one included already in (E4h. 
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A simple estimate with Mpi ~ lO^^GeV and M^, ~ 10~'^eV gives ~ lO^^mm. At 
the distance smaller than we will observe the four-dimensional world. Thus, this 
model predicts simultaneous modification of gravity at a distance of the order of 0.1 
millimeter and at around the Hubble distance! 

Let us now study the case of two sources placed at different positions y and y' in 
extra dimensions. The = 1 theory was studied in [^] (see also ||2^). In the present 
case the equation for the Green's function in the momentum four-space takes the 
form: 



(mi [(/ - a?) + ifAp'- SI 



The solution of Eq. 



+ M^i5(^)(y)/j G{p,y.y') 

= s^'\y-y'). 

can be written as follows: 



(34) 



G{p,y,y') 



1 



^ e-p\y-y'\ _ ^~\y~y'\\fi^+W^ 



4vr \y - y'\ 



M5 



M5 + M|ip2 
Using this expression one finds the potential: 




(35) 



y{r,y,y') 



G{t, X — x', y, y') dt 



dp 2p sin(pr) 



(27r)^ 



G(p,y,y') 



where G(p, y, y') is given by Eq. ( |5BD . All the expressions above should be un- 
derstood as an expansion in powers of [p'^/M^). The expression (^) reveals that 
the bulk interactions near the brane are strongly affected by the presence of the 
brane-induced terms. This is similar to what was observed in the A^ = 1 case in 
Ref. i. 

From Eqs. (^,^), it is clear that the potential between two sources placed at 
different positions in extra dimensions does not give rise to any ultra-local interac- 
tions 



5 Propagators on Thick Branes 

In this section we return to the consideration of branes with a finite width. We will 
check whether the results obtained for a delta function brane can be reproduced 
from a thick brane by the limiting procedure when A — 0. 
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The equation for the Green's function which is obtained from ( p4D takes the 
form: 

V Gip,y) + Ml,p'f\y) J d'^y' f\y') G{p,y') = g\y) , (36) 

where and are two locahzed functions with different thicknesses (we drop tilde 
sign over the momentum space Green's functions). V stands for an operator in 
dimensions which, as in the previous section, is approximated by the standard 
d'Alambertian plus one leading HD term: 

T> = |/ -V% + ^{p'- } . (37) 

Before we turn to the calculation of the Green's function let us make comments 
regarding the functions / and g. The former determines the localization of the 
breathing mode on the brane worldvolume while the latter is a measure of local- 
ization of a source under the consideration. In general these functions might be 
different. On the other hand, if we turn to the limit M^, — > 0, then Eq. (^) makes 
sense only if f ^ g in that limit. Thus / and g cannot be completely arbitrary func- 
tions, there should exist some correlation between them. This is a self-consistency 
requirement which will play an important role below. 

We look for a solution of Eq. (|36| ) in the following form: 



Gip.y) = Goip,y) + xiP.v) , (38) 
such that Go{p, y) obeys the equation 

VGo{p,y) = g\y), (39) 
and, hence, can be written as follows: 

GoiP, y) = j d^'y' g\y')D^{p, y-y'). (40) 

Here = D/M^~^^ , denotes the Green's function for the operator "D [D is defined 
in (^)). Note also that in the above expression for Go we dropped the zero-mode 
solution of the operator V which in this case is unphysical. 
The equation for the function x(p, y) takes the form: 

Vx{p,y) +Ml,p' [Y, + Y^] f{y) = , (41) 

where we introduced the following notations: 

Y, ^ Jd''y'f{y')Go{p,y'); Y, ^ J d'^y' f{y') x{p,y') ■ (42) 
Disregarding the unphysical solution of the homogeneous equation we obtain: 

Xip^y) = Ml, p' [Fi + Y2] I d'^y'fiy') D,{p,y - y') . (43) 
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In order to determine Y2 we integrate this equation with the weight Using the 
resulting expression we derive the following equation for x'- 

. ^ _ [Id^y'Piy')G,ip,y')] [J z f{z) D^p^y - z)] 

^^^'^^ 1 + M^y Jd^yd^y' Piy)D,{p,y-y') Piy') ' ^ ^ 

Therefore, the propagator in the general form reads as follows: 

G{p,y) = Go{p,y) - . J'J, ' Y{p,y) , (45) 

1 + Mpi 



where we have used the notations introduced in ( ^21) and 

Y,^j d'^y d^'y' f{y) D^p, y - y') f{y') , 

y (p, y) ^ j d^'y' f{y') D,{p, y - y') . (46) 

Below we will use these general expressions for the Green's function to derive 
the properties of the potential on the brane. 

Let us study the D = b case first. We parametrize the functions / and g as 
follows: 

f^y)=[^ forM<A, {A- forH<a, 

^ \ for |y| > A ; ^ \ for \y\> a . ^ ' 

First, for simplicity we assume that a < A. Then, inside the brane [y < A) we find 

_ 1 (. _ -pA sinh(pa) \ 
' 2AM3p2 1^ pa ) ' 

y ^ 1 (. _ g-pA sinh(pA)\ 
^ 2AM3p2 pA ; ' 

1 

2AM3p2 

At y = the first order correction to the propagator reads 

As we discussed above (in a paragraph after Eq. ( p6D the self-consistency of the 
equations for the Green's function (^) requires that A = a + 0(M*/M|i). This 
was used above and the leading term was kept. Thus, the correction due to the 
4D propagator appears at the scale A. However, there are other corrections due to 
the terms of order (9(M^,/Mpj) which distinguish a and A. These terms can give 
rise to the modification of the propagator at the momenta above M^. The latter 



Y^y) = , fl - e"^^ cosh(py)) . (4J 
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modification is expected to be in the theory anyway since the higher dimensional 
operators are suppressed by M^. 

As a representative example of higher dimensions we again concentrate on the 
codimension three case, = 3. For the functions / and g we choose the following 
interpolation: 

ny) = I /°'nf '^A ' 9\y) = I ^ /°',lf ^ " ' (50) 

■' \ for \y\> A ; ^ [0 for \y\> a . ^ ' 

As before we set a < A. Somewhat lengthy calculation gives the following result 
for the Green's function expanded in powers of the brane width A: 

1 / 3Q \ 

As in the A^ = 1 case we used the consistency relation A = a + 0{M^/M'^]). The 
corrections to the 4D behavior arise at the scale of order 1/A. However, this might 
be changed if we were to take into account a possible difference between A and a 
which is of order 0(M^,/M|j). The latter results in the corrections to the propagator 
which become significant at the scale M*. Hence, the propagators in D = 5 and 
D > Q are rather similar. 



6 The Crossover Scale 

In this section we will summarize the results on the crossover scale Tc- This quantity 
is determined by the 4D Planck mass Mpi and the bulk Planck mass M*. As we 
discussed before, in a 5D theory with a brane-induced Ricci term the expression 
takes the form: 

(D=5) _ . (52) 

While in the case of higher dimensions with zero-tension branes and the HD opera- 
tors suppressed by powers of we get: 



(D>6) ^ ^ (53) 



Thus, in higher dimensions (i.e., for -D > 6) the value of r^. is smaller. 

Let us now discuss what are the phenomenological bounds on Tc? The bounds 
could come from different sources. Let us start with the upper bound. There is 
a lower bound on the value of M* which comes from the data on sub-millimeter 
gravity measurements [jl| and the accelerator, astrophysical and cosmological data, 
that is M, > 10"^ eV 0]. Hence, the upper bounds on Tc are: 

^(D=5) < ^q59 . ^iD>6) < 1Q29 (54) 
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As we discussed in the introduction, the latter number coincides with the present 
day horizon size. Therefore, for D > Q the model with ~ 10~^ eV predicts 
simultaneous modification of laws of gravity at short distances around 0.1 mm and 
at large distances around the Hubble scale. 

Let us now turn to the lower bound on the crossover scale. This bound can come 
from the measurements of the Newton force at macroscopic distances, as well as 
from cosmological considerations. 

The distance at which Newton's law is known to hold exceeds somewhat the solar 
system size. Beyond this distance scale, because of the presence of dark matter, the 
Newton law can be modified if the simultaneous changes are made in the amount 
and distribution of dark matter so that these changes render a theory consistent 
with the data on the large scale structures. Therefore, from these arguments: 

r^"P > 10^^ cm . (55) 

However, the cosmological consideration can impose a stringer bound. Indeed, if 
the cosmological evolution changes at distances of order Tc, then we want to be 
of the order of the Hubble size at least. The cosmological solution of the 5D theory 
which changes the regime as the Hubble parameter becomes of order Tc was found 
by Deffayet in Ref. 0. Therefore, we have to impose r^^P > 10^^ mm. In this case, 
we get 10~^ eV < < 10'' eV in the 5D theory. However, in the D > 6 case we 
obtain a very stringent restriction, ~ 10~^ eV. Thus, the fundamental scale of 
gravity in this model is bound to be M^, ~ 10"'^ eV. 

Before we turn to the next section we would like to make some comments. Re- 
cently new interesting branches of cosmological solutions were found by Dick H,P|, 



and by Cordero and Vilenkin |]T0|. The creation of the "stealth" branes were dis- 
cussed in |]T0[. These solutions, from the perspective of a 4D braneworld observer, 
coincide with the well known solutions of pure 4D theory. On this branch, there is 
no difference between the time evolution of the universe in just a purely 4D theory 
from that in a higher dimensional theory discussed in the present paper. These solu- 
tions are found as follows [p|-[TO[|. Consider the Einstein equations with the induced 
terms included: 

= VsT^r" s^'^Ky) ■ (56) 

There can exist a solution which is fiat in the bulk but satisfies the equation on the 
brane, i.e.: 

TZab = 0, 

t2 f r> „ r>\ rj-iBrane 



The solution of the second equation in this system is exactly the one of a pure 4D 
theory. It is also a solution to the whole 5D system provided that the metric is flat 
in higher dimensions (see Refs. ||10[| ). 
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If this branch is reahzed, then there is no bound on Vc arising from cosmology. 
The only restriction on comes from the solar system measurements of Newton's 
law (0). This can be translated into the following upper bounds on M^,: 

10-3 eV < M, < 10^2 gv , for D = 5 ; 

IQ-^ eV < < 10^ eV , for D > 6 . (58) 

As before, there is a bigger range for in the 5D theory. 



7 Comments on Nonzero Tension Branes 

In the previous sections we dealt with branes (or system of branes) which were 
placed at orbifold fixed points and for which the tension effects were removed by 
fine tuning. Here we would like to study the influence of a nonzero brane tension on 
the expression for r^. We start with an infinite (4 + iV) dimensional space. Consider 
a certain brane in this space which has an intrinsic tension T. The brane distorts 
the ambient space around it. For a BPS D3-brane the solutions is regular outside 
of the core and has well defined horizon. However, for non-BPS branes which are 
phenomenologically relevant, the solutions give rise to singularities. The singularity 
at the core is an artifact of the delta-function approximation is expected to be 
removed if the core of the brane is smoothed (see also [0). We will adopt 
below this philosophy. 

For a brane world which preserves 4D Lorentz invariance the interval takes the 
form: 

ds^ = A\y) r^^.dx^dx'' - B\y) dy' -C\y)yUnl_, , (59) 

where the functions A, B and C depend on the brane tension and M^,. For certain 
cases A^B and C are known exactly ( see, e.g., Refs. |30|,|32|,|33|). In the N = 2 
case a positive-tension brane creates a deficit angle and produces no other nontrivial 
warp-factors like A, B and C. The effects of this on the model with brane-induced 
Einstein-Hilbert term were studied in Ref. where it was concluded that the 
qualitative behavior found in (and outlined above) does not change. 

Below we will concentrate on the case > 3. We use the solutions for non-BPS 
branes obtained in Ref. [pO|] : 



_1 I N-l 1 ( iV-3 / N-l \ 1 ( 1 . I N-l \ 

where 

/-I + (^) • (61) 

Here the gravitational radius of a brane in the transverse direction yg is determined 
by the brane tension T and can be written as follows 
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In order to determine the properties of the interactions on the brane one needs to 
calculate a graviton propagator on the background (|60D. The exact calculation is 
cumbersome (for certain calculations along these directions see Ref. [0). However, 
one can do a qualitative analysis of the equation for perturbations. The analysis is 
similar to that performed in Ref. ||^ and will not be repeated here. We just mention 
the basic steps. 

The solutions for A, B and C in (|60D are singular near y = 0. However, in this do- 
main the classical theory breaks down and the expressions ( |60D should be smoothed 
out. In the previous section we accounted for this UV regularization by introducing 
higher dimensional operators. This is perhaps the most consistent and reliable way of 
dealing with the singularities. However, when a nonzero brane tension is taken into 
account, the higher dimensional operators complicate the equations even further so 
that they are not treatable analytically. Under these circumstances, it becomes more 
reasonable to introduce other UV regularization, e.g., one could use the finite width 
of the brane core A. Then, we can estimates the critical momentum/ distance for 
which the 4D induced term becomes dominant over the bulk terms in the background 
(|60|). Note that this approach has a rather limited applicability since it is justified 
only when T ~ M^. This consideration leads to the following region where the ef- 
fects of the 4D induced terms dominate: r < rf-^''*) ~ Mpi/^Jr + +2. 
Here the 4D Planck mass is determined by the product of the induced scale and the 
value of the function A taken in the core M|i = M{^^A^(A). Although, the value 
of A^{A) is not known, one might expect that it is not far away from the value 
of this function at the point where the gravity approximation is reliable. Thus, it 
adds/subtracts an order of magnitude from Mind- Note also that since for a phe- 
nomenologically viable model M^, ~ 10^^ eV, then the expression for r^^"*^*-* is 
applicable only for T < (10~^ eV)'^. 

Let us now turn to the constructions in which our brane is represented by a 
stuck of branes and anti-branes (or orientifold planes) which can be placed at some 
orbifold fixed points. In this case one can study consistently the limit T — > 0. 
The resulting expression for r^^^'^*) coincides with that obtained in Ref. 0. The 
crossover scale depends in general on the number of dimensions. On the other hand, 
in the previous sections we studied a theory with high derivatives and obtained 
the crossover scale which was independent of (^3|). The question is how these 
two results are reconciled. The reconciliation is based on the following arguments. 
When we take into account the higher dimensional operators which are suppressed 
by M,,, this means that UV gravitational resolution of the brane width is at the 
scale A ~ Substituting the latter expression in the expression for r^""^'^*^ 

with T = we recover Eq. (p3[). Thus, the phenomenological constraints on the 
bulk gravitational scale are the same as in the previous section. 

Finally, the following scenario is also possible: suppose the bulk gravitational 
scale is high, let us say in a TeV region or so. Let us call this scale Msuik ^ TeV. 
It is conceivable that some scalar field dynamics suppresses the coefficient in front 
of the Ricci scalar in the bulk. Then, Mt+^ = /((<l>))M^+i^, where /(($)) is a 



20 



suppression factor which depends on the vacuum expectation value of a scalar. On 
the other hand, the higher derivative terms in the bulk are suppressed by the scale 
-^Buik ^ M*- Let us consider a zero-tension brane at the orbifold fixed point. Then, 
in the expression for r^"^~^^*) we should take A ~ I/MbuiIc- Using the hierarchy 
-/^Buik ^ we can obtain an acceptable values of Tc even for ~ TeV. For 
instance if Meuik ~ ^pi we get the value of which is consistent with (^5D for 
AT > 4. 

8 Conclusions 

Summarizing, we discussed a brane-world model with two or more infinite-volume 
extra dimensions. This model predicts the simultaneous modification of the law of 
gravity at around 0.1 mm and at the Hubble scale. Therefore, the fundamental scale 
of gravity is bound to be 10"'^ eV in this model. These predictions can be tested by 
table-top gravitational measurements and by astrophysical experiments. 

The UV singularities which arise in this model can be treated consistently. We 
used both, the regularization by a finite brane width, and by the higher-dimensional 
operators which are suppressed by the fundamental scale of gravity. We calculated 
the corresponding two-point Green's functions in the regularized theory and showed 
that there is an intimate relation between the UV and IR scales where gravity 
changes its regime. 
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